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I. INTRODUCTION 

Since the original work by Fisher and de Gennes [Tj, 
it has been well accepted that in the vicinity of a critical 
point the developing long range fluctuations give rise to 
a Casimir force in a system with a film-like geometry. As 
a consequence of finite size scaling [5], the generic form 
of the force is 



Casimir 

ksTc 



(1.1) 



where Tq is the critical temperature, L is the thickness 
of the film, d is the bulk dimensionality of the system, 
^ is the correlation length which diverges at the critical 
point in the bulk, and i9 is a universal function of its one 
argument. 

The notion of universality requires some qualification 
in the context of Casimir forces. Generally, for a sys- 
tem with a critical point, universality implies dependence 
only on the universality class to which that system be- 
longs. In the case of an Ising model, that universality 
class is characterized as 0(1), corresponding to the dis- 
crete symmetry that is broken in the ordered phase. Sys- 
tems belonging to this class include uniaxial ferromag- 
nets, simple liquid- vapor systems and systems exhibiting 
de-mixing [3 . 

The universality class, however, is not the only char- 
acteristic of the system that dictates the behavior of the 
function i?(a;); one must also consider the boundary con- 
ditions at the two extremities of the film. This effect 
is dramatically exhibited in the results of measurements 
performed in Ref. |3J H]. Those experiments monitor 
the thickness of a film of "^He as the temperature passes 
through the lambda point at which bulk Helium acquires 
superfluid characteristics. It was found that the thickness 
of the film, which directly tracks the Casimir force ap- 
pearing in the vicinity of this critical point, changes sub- 
stantially when the temperature is close to and slightly 
below T\. Quite interestingly, deep in the superfluid 



phase, the film partially recovers its thickness but still 
remains thinner than in the normal phase. 

The thinning of the film well below the transition point 
was successfully explained in Ref [B], where it is shown 
that the confinement of bulk Goldstone modes combined 
with the effect of finite film thickness on surface fluc- 
tuations gives rise to a Casimir force that reduces the 
thickness of the film in the superfluid phase as compared 
to that of the normal film. However, the much more sub- 
stantial thinning of the film right below the critical point 
remained unexplained for many years. 

Earlier predictions of the Casimir force at the vicinity 
of A transition based on calculations emphasizing the cor- 
rect universality class, 0(2) for the superfluid transition 
in ^He [ZIIH], differ markedly from the experimentally de- 
termined force, which is greater in magnitude by a factor 
of about thirty than those theoretical results. 

In Refs. [9l-fTT| it was shown that the discrepancy arises 
from the boundary conditions assumed in the previous 
calculations. In particular, prior field theoretical calcu- 
lations [21 [12] as well as the simulations in [2 assume 
periodic boundary conditions, while physical reality re- 
quires the superfluid order parameter to vanish immedi- 
ately outside of the Helium film, thus mandating Dirich- 
let boundary conditions. It has been known for a long 
time that the boundary conditions operative in a sys- 
tem play an important role in determining the behavior 
of Casimir forces [131 [H]. However, there is no a priori 
reason to assume that the strength of Casimir force with 
Dirichlet boundary conditions will be significantly larger 
than that for the case of periodic boundary conditions. 
In fact, the absence of this expectation undoubtedly hin- 
dered our understanding of the dip in Helium films for 
many years. 

The strong effect of boundary conditions on the 
Casimir force in the particular instance of an 0(2) system 
in three dimensions leads one to ask if there are any gen- 
eral inferences one can arrive at regarding the interplay 
of boundary conditions, dimensionality, and universality 
class on the Casimir force acting on a system at and 
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near the critical point. With an eye toward developing 
a greater depth of understanding of effects of boundary 
conditions and dimensionality, we have looked closely at 
one system for which exact results can be obtained: an 
Ising model in one and two dimensions. With regard to 
the one-dimensional Ising model relevant results follow 
almost immediately from a simple transfer matrix anal- 
ysis. In the case of the two dimensional Ising model, we 
have been able to perform an analysis of that model that 
yields the critical Casimir force for four different bound- 
ary conditions: periodic, antiperiodic, free (the appro- 
priate surrogate for Dirichlet boundary conditions in the 
Ising model) and fixed-spin. The results of this analysis 
are consistent with those we have extracted from various 
expressions already in the literature. 

A key conclusion that we draw is that boundary con- 
ditions alone do not determine the relative amplitudes of 
critical Casimir forces. Indeed, we find that as dimen- 
sionality is reduced, periodic (and antiperiodic) bound- 
ary conditions give rise to a critical Casimir force that 
is increasingly stronger than the force resulting from free 
boundary conditions. In fact, the critical Casimir force 
for a one dimensional Ising model with free boundary 
conditions is effectively non-existent. 

The paper is organized as follows. In section [ll] we 
calculate the Casimir force for the one-dimensional Ising 
model as a prelude to the discussion of the two dimen- 
sional version. We obtain the Casimir forces with free, 
fixed, periodic and anti-periodic boundary conditions. 
We observe significant differences in the forces in these 
four cases. In section \LU\ we present our results for the 
two dimensional Ising model and compare them with the 
one dimensional results. The forces are derived in two 
different ways. First, we show how they follow from ex- 
pressions already in the literature. We also demonstrate 
how they arise from an unified analysis of the two di- 
mensional Ising model making use of an approach along 
the lines of the method developed by Schultz, Mattis and 
Lieb [15 . In section |vj we discuss our findings, their im- 
plications, and summarize our conclusions. For clarity 
of exposition, details of our calculations are relegated to 
appendices. 



II. ONE DIMENSION 



A. periodic boundary conditions 



The spin variables Si can assume the values ±1. Given 
the transfer matrix 



(2.2) 



where /3 = l/ksT, a spin state at the left boundary of 
the A^-spin array of the form {l\ and a spin state \r) at 
the right boundary of that array, the partition function 
is equal to 



(2.3) 



The quantity on the right hand side of (2.3) is most con 



veniently evaluated making use of the eigenvector decom- 
position of Ti 



Ti = |e)2cosh(/3J)(e| |o)2sinh(/JJ)(o| (2.4) 



where 



1 

1 

71 



(2.5) 
(2.6) 



When periodic boundary conditions hold, we calculate 
the partition function by evaluating the trace of Tf over 
a complete set of bounding spin states. Given (2.5) and 
(2.6) we have immediately 



odic(Af) 



cosh^(/3J) +sinh™(/3J) , (2.7) 



N, 



Therefore, the free energy of this system is 



= - A^fcsT In (2 cosh(^ J)) - fcsT In [l + tanh(;9 J)^] 

(2.8; 



In the thermodynamic limit the second term of ( |2.8[ ) be- 
comes negligible, so the bulk free energy per unit length 
(alternatively per spin) is 



fs = -ksT \n2- kBT\ncosh{l3J) 



(2.9) 



The correlation function, C(Zi,^2), and the associated 
correlation length, ^, in the one-dimensional Ising model 
are defined as follows [11] 



tanh(^J)'2-'i 
exp[-(/2-/i)/C] 



(2.10) 



which leads to 



The partition function of the one dimensional Ising 
model is well known to be derivable from the two-by-two 
transfer matrix of that model. The Ising model possesses 
a global Z2 symmetry and is described by the following 
Hamiltonian 



(2.1) 



where J > is the spin-spin coupling constant, and the 
sum is over nearest neighbor bonds z, j on the lattice. 



C = -l/(lntanh/3J) 



(2.11) 



Taking the derivative of the free energy ( 2.8 1 with respect 



to N and subtracting the bulk free energy per unit length, 
we are left with 



ON 



+ fB ^ 



^ l + e^/« 
knT (N e-^/« 



A^ 



^ 1 



77'CaS 

Id per 



(2.12) 
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The Casimir force has precisely has the desired form 
(1.1). That is, for the one-dimensional Ising model 



^id per (a;) 



1 



(2.13) 



Note that the factor fc^T plays an essential role in our 
result for the Casimir force. To the extent that it can 
be said to exist, the critical point of the one-dimensional 
Ising model is at T = 0. Unless we normalize the criti- 
cal Casimir force in terms of this temperature-dependent 
factor, we are forced to conclude that it vanishes. 



B. Antiperiodic boundary conditions 




^Id per{x) 



FIG. 1: The universal Casimir force function -did pcr(a;), 
shown as a dashed curve and the function -did antipcr(a;), 
shown as a solid curve. 



The second boundary condition we study here is an- 
tiperiodic. We enforce this boundary condition by intro- 
ducing an array of antiferromagnetic bonds in the system. 
In this case, there are A'^ — 1 two-by-two transfer matrices 
of the form (2.2) and one transfer matrix of the form 



'0J 



(2.14) 



The eigenvector decomposition of the transfer matrix T2 
is 

T2 = |e)2cosh(/3J)(e| - |o)2 sinh(^J)(o| (2.15) 

The partition function of this syster n is the tr ace o f the 
uant 
have 



quantity Tf • T2. Making use of (|Z5|) and (|2J|), we 



-'antipcr 



Tr (Tf • T2) 

(2 cosh pjf - (2 sinh /3 J)^ (2.16) 



This means that the free energy is given by 

J" = -kBTlnZ 

= -NksT In (2 cosh /3 J) - fcsT In [l - (tanh /3 J)^] 



NfB-kbT\n 



(2.17) 



The quantity fs in the last line of (2.17) is the "bulk" free 



energy of the one-dimensional Ising model. The quan- 
tity ^ = — l/(lntanh f3J) is the correlation length of that 
model. We then take the negative derivative of the free 
energy with respect to L, add in /b, and end up with 



iTiCas 
Id antipcr 



fE 



e l-e-^/« 
ksT N e-^/« 

keT 
N 



^id antiper (iV/0 (2.18) 



Fiture [T] displays universal functions -did pcr{x) and 
"^id antipcr (a;)- Two distinct features of this scaling func- 
tion distinguish it from the case of periodic boundary 
conditions; it is positive, indicating a repulsive Casimir 
force, and it does not go to zero at the zero temperature 
critical point. 



C. Free and Fixed Boundary Conditions 

In the case of free boundary conditions, the partition 
function of the A^-layer system is given by 



Zfrec = 2(e|T^ 



2(2cosh(/3J)) 



N 



(2.19) 



Taking the derivative with respect to N of the log of 
the above result and subtracting the bulk contribution 
to that derivative, we find perfect cancellation. There is 
no residual force. 

When the spins are fixed at the boundary, the state 
between which the N^^ power of the transfer matrix is 
sandwiched is 



I/) 



(2.20) 



and the partition function is given by 

-Zfixcd 

= (/|Tf I/) 

= ^ ((2cosh(/3J))^ + (2sinh(/3J)) 



(2.21) 



Given the near identity between the final line of (2.21) 



and what was obtained for the partition function of the 
system with periodic boundary conditions in (2.7), we 



find immediately that the critical Casimir force induced 
by fixed boundary conditions matches the corresponding 
force in the case of periodic boundary conditions. 
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III. TWO DIMENSIONS 

In this section, we consider an Ising model on a two 
dimensional square lattice consisting of M rows and N 
columns. The Hamiltonian of this model is the same as 
given in Eq. (2.1 ) but the sum is over a two dimensional 
lattice. We will begin by focusing on periodic and an- 
tiperiodic boundary conditions, first demonstrating how 
the critical Casimir force follows from expressions for the 
partition function already in the literature. Then, we 
perform an analysis of the two dimensional Ising model 
based on the fermion variable approach introduced by 
Schultz, Mattis and Lieb that leads to precisely the same 
results. Finally, we make use of this same approach to 
rederive existing expressions for the critical Casimir force 
in the case of free and fixed boundary conditions. As has 
already been noted, the Casimir forces in the two di- 
mensional Ising model for free and fixed-spin boundary 
conditions are connected by a duality relationship. 



A. Periodic boundary conditions 

Our starting point in the calculation of periodic bound- 
ary conditions is the formula utilized by Ferdinand and 
Fisher [T7] for the partition function of a M x TV square 
Ising lattice wrapped on a torus. 



exp 



'N-l 



+ exp 



^ In cosh ^Af 72^+1 

N-l ^ 

In cosh - Mj2r 



exp 



'N-l 



' exp 



^ Insinh ^My2r+i 
Insinh ^Mj2r 



r=l 
'N-l 



In (3.1) 



li = In 



ci + ^/cf - 1 



(3.2) 



where 



Q = cosh 2i4r coth 2 AT — cos ^Tr/n 
= cosh2A:coth2A: - + ^) 
= c{zi) 



with 



and 



K = 



J 
kef 



zi ^ e 



Hit / n 



(3.3) 



(3.4) 



(3.5) 



Given that we are interested in the limit N ^ oo, 
corresponding to the thermodynamic limit for the width 
of the of the two dimensional strip with finite thickness 
A/, we can replace the sums in the exponents of (3.1) 



by integrals. This means that the partition function in 
question reduces to 



2 exp 



N 
2^ 



2ir 



2 exp 



N 
2^ 



In sinh 

27r 



M 



lie'*) 







In cosh ( Y7(e*'^) 



(3.6) 



(3.1) which can be rewritten as follows: 



2 exp 



'N 
2^ 



In cosh ( —-f{e"^) 







X 1 -I- exp 



N 
2^ 



2 exp 



' N 
2^ 

2ti 



2 

Insinh I — 7(e*^; 







In cosh ( y7(e''^) ) dcj) 



In cosh ( ^7(6*"^) 







1 -I- exp 



N 
2^ 



2ir 



Intanh ( y7(e'"^) 



(3.7) 



Because tanha: < 1 for any positive x, the function 
Intanh a; will be negative. Thus the integral on the last 
line of (3.7) produces a negative number, which is multi- 



plied by N, the transverse extent of the strip, the analog 
of the area of a three dimensional film. We are interested 
in the limit n — cx), so the second term in curly brackets 



on the last line of (3.7) will be negligible compared to 



unity, and the sinh contribution to the partition function 
in (3.6) can be discarded without error. 



Taking the log and then the derivative with respect to 
m, we end up with the following integral 



a J" ksT 



DM TT 



7(e**)tanhAf7(e''^)d0 (3.8) 



Making use of the fact that limj,_>.oo tanhx = 1, we cal- 
culate the force for the bulk system which we subtract 
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from (3.81 to find the following Casimir force 

l92dpor 



-[ j{e''^) tanh M-f{e"^)d(j> - [ j{e"^)d(j> 

Jo Jo 

(3.9) 



B. The scaling limit 

To obtain the scaling form of the Casimir force, we 
start with looking at the behavior of the function 7(e"^) 
in the vicinity of the critical point at which 



cosh 2X coth 2is: = 2. 



(3.10) 



The solution to the equation (3.10) is given by 

sinh^Kc = 1 

Kc = ^ln(l + ^(2)) 0.44068... (3.11) 

as given in Ref . jl7] . We now introduce the reduced tem- 
perature, r: 



T = - sgn ( i cosh 2K coth 2i^ - 1 



X A / - cosh 2K coth 2K - 1 



(3.12) 



Then for small r and 0, c(e^'^) = 1 + 2x^ + 4-, and for the 
corresponding expansion of the function 7(0(6**^)), we are 
left with 



7(c(e'^)) 



In 



1 + 2v/t2 + </'2/4 



2v/r2 + 4,y4 



(3.13) 



If we set (p = 2uj and replace lu by fl/AI and r by x/M, 
we find the following expression for the Casimir force t 



1 

7rM2 



dn^/x^ + [tanh y/x'^ + - 1 



(3.14) 



The universal function •d{x) that emerges from (3.14) is 



(x)^ - dVL^Jx'^ + tanh v/^^+TF - 1 



(3.15) 

This result is displayed as the short dashed curve in 
Fig. [I 



C. The two dimensional Ising model with 
antiperiodic boundary conditions 

The exact partition function for a two dimensional 
Ising model on M x N square lattice with anti-periodic- 



antiperiodic boundary conditions is 



exp 



N-l 



■ exp 



^ lncosh^M72r+i 

L 

N-l ^ 

lncosh-M72r 



r=l 
'N-l 



r=l 



exp 



exp 



JV-l 



^ lnsinh^M72r+i 

r=l 

In sinh - M72. 



(3.16) 



The derivation of the above partition function is given in 
Appendix [X] and is based on the expression calculated by 
Wu and Hu in Ref. [T5] . Note that the last term in ( 3.16 ) 
changes sign at the critical temperature; see Appendix |A| 
for more detail. As previously, quantity M is the width 
of the strip, while N is its extent. Again, we will take 
the limit TV — 00 before allowing M to become large. 

The difference between periodic and antiperiodic 
boundary conditions lies in the sign relationships be- 
tween the terms in Eqs. (3.1) and (3.161. Consider the 



two terms with the log of the hyperbolic cosine in Eq. 
(3.16). They appear with opposite signs, and as we will 



see, this leads to near-cancellation of their contributions 
to the partition function. On the other hand, the other 
two terms in (3.16) may or may not add up. In order 



to assess the full partition function, we need to look at 
what remains after all near cancellations are taken into 
account. 

Appendix[B]contains an analysis of the contributions of 
the various terms in (3.16) to the partition function. The 



end-result is that the hyperbolic sine terms now domi- 
nate. We take a derivative with respect to M, the film 
thickness, to obtain the Casimir force. Performing the 
same scaling analysis as in Section p^II A[ we find for the 
Casimir force in the case of antiperiodic boundary con- 
ditions. 



1 



7rM2 



dn^/^ 



02 



coth + fi2 _ I 



(3.17) 



For comparison, see (3.14) for the Casimir force in the 
case of periodic boundary conditions. Figure [2] shows the 
Casimir forces for periodic, antiperiodic and free bound- 
ary conditions. Note that when antiperiodic boundary 
conditions apply, the critical Casimir force is repulsive, 
leading to film thickening, as was the case for the one 
dimensional Ising model. 



D. Free boundary conditions 

To obtain the Casimir force for a two dimensional Ising 
model in the case of free boundary conditions, we use the 
expressions used in Li, et. al. |19j . After some manipula- 
tions, we find the asymptotic scaling form of the critical 
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portion of the free energy per unit area is 



-^singular 

1 /"^ duj 



2\n{l -Zsf 



+ ln 
1 



N 



+1 



1 - 2^/T^ 



N 



(3.18) 



Here Zs = taiih{/3Js) with J5 the in-layer couphng 
on the surface and r is the reduced temperature. In a 
change of notation, we now denote the thickness of the 
film by N. For simplicity, we assume that bulk couplings 
in both X and y directions are equal. Note that the re- 
duced temperature in the above equation is equal to 1 /N 
of the reduced temperature defined in Li et. al. [H]. To 
find the Casimir force, we take the negative of the deriva- 
tive of the free energy above with respect to N. Since N 
appears only in the exponent, we need to consider the 
following derivative 



_d_ 
dN 



In (1 + 



N 



N 



(3.19) 



which helps us introduce the relevant scaling variables. 
If we replace the reduced temperature, r, by x/N and 
the integration variable, ui, by fl/N, the right hand side 
of (3.19) in the limit of large N becomes 



N 



(3.20) 



and thus the negative derivative of the free energy (Eq. 
Slsl) is 



dR 



singular 

~dN 



dn^/x^ + Q^x 



1 



^2Vx^+n^ + ( 1 - 



Vx'^+n^ 



^/x■^+n^ 



'2Vx'^+n^ 



-2^x'^+Q? 



(3.21) 



Now we need to subtract from the above expression the 
"background" casimir forces which is obtained when the 
thickness of the film is infinite. As TV — >■ 00, the free 



energy given in (3.181 becomes dominated by the term 



going as (l-H2r/-\/T^ -l-w^)^, and thus the corresponding 
Casimir force per unit area is 



1 



djjj 



lnfl-F2v/^ 



(3.22) 



Inserting the scaling variables x /N = r and uj — VL /N as 
before, we obtain 



-4t^ I dfly/x'^ + (3.23) 



Subtracting this from (3.21), we are now allowed to ex- 



tend the upper limit of integration to 00, and thus we 
end up with the final expression for the Casmir force per 
unit area 



ttN^ 



,2Vx''+n^ 



-2Vx^+n2 



- 1 



(3.24) 



which is plotted in Fig. [2] (long dashed line) and as ex- 
pected is negative giving rise to an attractive Casimir 
force. It is interesting to note that the strength of the 
Casimir force in the case of free boundary condition is 
zero in one dimension but in two dimensions is compa- 
rable in amplitude to the Casimir force with periodic 
boundary conditions. We also anticipate the results of 
the next section by plotting the Casimir force for fixed 
spin boundary conditions, previously obtained by Evans 
and Stecki [20], which is related to the Casimir force for 
free boundary conditions by duality. 



IV. ALTERNATIVE APPROACH TO THE TWO 
DIMENSIONAL ISING MODEL 

It is also possible to obtain all the expressions above by 
a set of direct calculations of the partition function of the 
two dimensional Ising model in zero magnetic field. We 
make use of the method introduced by Schultz, Mattis 
and Lieb [T^ , which will hereinafter be referred to as the 
LSM method. The version of this model we will address is 
wrapped on a cylinder with free boundary conditions on 
the ends. Circumferentially, we will allow for either pe- 



7 



or4 

/ 0.2 






K/ 



X 



FIG. 2: The scaling contribution in the two dimensional 
Ising model to the Casimir force for antiperiodic boundary 
conditions (solid curve), periodic boundary conditions (short 
dashed curve) and free boundary conditions (long dashed 
curve) and fixed spin boundary conditions (alternatively long 
and short dashed curve). 



riodic or antiperiodic boundary conditions. In this way, 
we encompass all the boundary conditions that were con- 
sidered above. Additionally, we will demonstrate by a 
duality transformation that the critical Casimir force for 
free boundary conditions is simply related to the Casimir 
force for fixed boundary conditions obtained previously 
by Evans and Stecki [20] . 

Denoting the partition function for free boundary con- 
ditions by Zf, our first step is to write Zf in terms 
of a transfer matrix along the axis of a cylinder with 
circumference M and length N. In the representation 
with two dimensional Pauli spin matrix crj, diagonal, 
J = 1, . . . , M, the transfer operator within a row is given 

by 



v. 



exp 



M 

i=i 



(4.1) 



and between rows by 

Vi = exp 



M 



(4.2) 



In (4.1 1 and ( |4.2[ ), if„ with n = 1, 2 are, respectively, the 
vertical and horizontal couplings, both positive for the 
ferromagnetic case, and K* is the dual of K^, i.e. 



tanhiiT* = exp [-2Kn] 



(4.3) 



Writing Vi in the form (4.2 1 requires that we accompany 
it by a factor {2smh{2Kiy)^^^^ whenever it appears in 
an expression generating the partition function. In this 
article we make use of the following representations of 



the Pauli spin matrices ax and cr^ 

1 

-1 

-i 

1 



-1 
-1 



(4.4) 

(4.5) 

(4.6) 
(4.7) 



The cylinder is assumed to couple to itself, so the ac- 
tual depiction of the system ought to be as a torus, as 
shown in Fig. [3] Free boundary conditions are enforced 




FIG. 3: The M x N lattice with M = 9 and TV = 7, shown 
as a torus to indicate the cyclical nature of the vertical bond 
arrangement. The vertical bonds between the seven rows are 
shown dashed. 



by changing the strength of one array of vertical bonds to 
Kq and then by taking the limit Kq 0. The partition 
function Zf, is then given by 



(2sinh(2i^i))*^(^-i)/2 ijj^ (2smh2Ko)^'/^ 



xTr I (V2Vi)^-iV2exp 



M 



E < 



m— 1 



(4.8) 



Given the fact that Kq — > cx) as Kq 0, we find 

M 



lim e '^'^'^0 exp 



E < 



m— 1 



|0)(0| (4.9) 



where |0) is defined by 

a||0)=-|0), j = l, 



(4.10) 
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Furthermore, one can readily verify that as Kq — )■ 0, 
(2sinh2Xo) exp(2_firg) — )• 4. Cychcity of the trace then 
reduces (14. 81) to 



(2sinh2i^i)*^(^-i)/22^^ 
x(0|(V2Vi)^-iV2|0) 



(4.11) 

begin- 
ning and end of the product of operators in (|4.11| and 



X /2 1/2 

The final step is to insert Vj^ = 1 at the begin- 



(4.12) 



then use vY^lQ) = cxp{-MKl /2)\0) . This yields 

Zf = (2sinh2Xi)'''''^^'^/'2*^e-^'^^i 
x(0|(V')^|0) 

where the operator V is given by 

V = Vi/^zVy^ (4.13) 

Given that the operator V' is self adjoint, a natural 
step is to make use of its eigenvalue decomposition as de- 
termined by Schultz, Mattis and Lieb [TS]. Before doing 
this, we will show how the strip with "hard" walls ad- 
mits of a similar representation. We consider the same 
toroidal arrangement. This time, we alter the interac- 
tions between the spins in one of the M spin circles. 
That is, we replace the interaction strengths in one of 
the V2 terms by Kq where we now take Kq to infinity. 
In addition, we add a term to the Boltzmann constant 
that keeps the partition function finite in the Kq 00 



limit. Then, (4.8 1 is replaced by 

= (2sinh2i^i)*^^/2 ijj^ Tr I (ViVz)"^"' V 



exp 



Given that 



lim exp 



where 



M 



^0 ^ i'^m'^m+l - 1) 



Kr 



m—1 



M 
m— 1 



1) 



(4.14) 



= !+)(+! + |->(-| 

(4.15) 



±l±) 



m 



.M 



(4.16) 



Given that the |±) are also eigenvectors of V2, we are 
led to another version of (4.14): 



(2sinh2i^i)*^^/2e-*^^''{(-t-|V 



where 



V 



(4.17) 



(4.18) 



Finally, it is straightforward to verify that for periodic 
(or toroidal) boundary conditions 



Zt = (2sinh2i^i)'*^^/2r[,^.v 



N 



(4.19) 



It is instructive to consider what can be deduced by du- 
ality. We can avoid having to consider the effect on eigen- 
vectors by considering (4.8) and (4.14) and examine the 
pair of operator products before taking the trace or the 
limit. A unitary transformation of the Fermi operators 
can be constructed [21] which interchanges 
and cr cr J f7j_|_]^; there is a technical problem involving 
the boundary, which has been overcome elsewhere |21| . 
This converts (4.8) into (4.15), but with Kq replaced by 
Kq. The limit goes through and (4.17) is identical up to 
numerical factors, provided K^ and K2 are interchanged. 
Thus, if we have evaluated the right hand side of one or 
other of these two equations, the other follows directly. 
Duality in the usual form is discussed in the appendix. 
What we are using here is an operator version of it as 
an automorphism of the algebra of Fermi operators which 
appears to have been known to Onsager, although this 
work was not published, as far as we know. 



A. The LSM approach 



In the LSM method Fermi operators /„, fl are in- 
troduced by means of the Jordan- Wigner transformation 



fk = Pj-i<yJ, 2 < j < M and /i = CTf (4.20) 
where a~ is the spin lowering operator, and Pj is defined 

by 

j 

^j=n(-'^i)' (4-21) 

1=1 

It is easy to check that the fj and /j are, indeed, Fermi 
operators, and that 



Vi = exp 



M 



-KiY,iv]fj~^) 



(4.22) 



The operator V2 has a more subtle form: 

C A/-1 



V2 = exp 



-PmUI- fM){fl + h) 



One can verify that 



[Vi,2,Pm] =0 



(4.23) 



(4.24) 
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so that, for instance, V and Pm can be diagonalized 
simultaneously. The operator Pm has eigenvalues ±1. 
Thus, we can write 



where 



and 



V' = P+V'(+)+P_V'(-) (4.25) 
V'(±) = Vi/^V2(±)Vi/^ (4.26) 



{M-l 



V2(±) = exp 

T(/I,"/m)(/i^ + /i) 



(4.27) 



The operators P± are projectors onto the eigenstates of 
Pm- 

Since Pm\0} = |0), it follows that 
Zf = (2sinh(2ifi))^-^(^-i)/22Mg-MKr^O|(y/(^))jV|o^ 

(4.28) 

It should be noted that if a line of reversed K2 bonds 
were to be inserted on every tier of the cylindrical lattice 
between site M and site 1, then the resultant partition 
function, which we denote Z'^, is given by 

Z] = (2sinh(2ifi))*^(^-i)/22*^e-*^^i(0|(V'(-))^|0) 

(4.29) 

The next step is to introduce running wave vectors as 
follows: 



M 



(4.30) 



Both periodic (for the minus sign) and antiperiodic (for 
the plus sign) wave numbers are required to reduce (4.27) 
to the form 



V2(±) = exp 

iF^i~k)+F{k)) 



K2 e^HFHk) - F{-k)) 

kenMi±) 



(4.31) 



where the sums over k are to be taken in the range 

nM{±) = {k : exp[iMk] = =f1, -tt < k < vr} (4.32) 



Equation (4.31) can be analyzed further by grouping 



equal and opposite k values together, giving 



V2(±) - exp 



2K2 {t" (k) cos k + T^ik) sink} 

(4.33) 



0<fc<7r 



where 



and 



'{k) = {F\k)F{k) + F''{-k)Fi-k) - 1) (4.34) 



ry{k) = -i{F\-k)F''{k) - F{k)F{~k)} (4.35) 



Terms in ( 4.33[ ) for different values of the summand com- 
mute. In the paired subspaces, the t"(/c), a — x,y, z are 
indeed spin 1/2 operators. Thus, we can write V2(±) as 
a product 



V2(±)= Yi exp [2K2 {T^{k) cos k + Ty{k) sink}] 

(4.36) 



0<fc<7r 

Similarly, for the operator Vi 

Vi(±)= n exp[-2i^*r-(fc)] (4.37) 



0<fe<7r 

Thus, V'(±) will also factorize: 



v'(±)- n v'(fc) 



(4.38) 



0<fc<7r 



where the k values are determined according to (4.32 1 



The behavior at fc = and k = n must be specified 
carefully: 



V (0) = exp [iK2 - Kl) (2Ft (0) F (0) - l)] (4.39) 



and 



V (tt) = exp [- iK2 + Kl) (2F^ (vr) F (tt) - l)] (4.40) 

Then, the eigenstates of these operators are easily con- 
structed. For the remaining Y'{k), note first that 



V'(fc)F^(±A:)|0) = F^(±/c)|0) 



(4.41) 



In the subspace with vectors |0) and F^{—k)F'^ {k)\0) , the 
operators Ty''^{k) introduced in (4.34) and ( |4.35[ ), along 
with the operator 



r^(fc) = F'^{~k)F\k) + F{k)F{-k) 



(4.42) 



behave like Pauli spin operators, as was observed in |15) . 
This shows that for fc 7^ 0, tt, 

\'{k) 

= cosh(7(fc)) — sinh(7(fc)) 

X (cosh((5'(fc))T^(/c) - sin((5'(fc))T^(fc)) (4.43) 

where the function "/(k) is determined by 

cosh(7(A:)) = cosh(2i4:*) cosh(2i<:2) - 

- sinh(2ifi*) sinh(2X2) cos(fc) 

(4.44) 
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where ^{k) is also non-negative for real argument. The 
angle S'{k) is defined by 

cosh(2if2) = cosh(2X*)cosh(7(fc)) 

- sinh(2i^i ) sinh(7(fc)) cos((5'(fc)) 

(4.45) 



and 



sin((5'(fc)) sinh(7(fc)) — sin fc sinh(2i4r2) 



(4.46) 



Notice that the term multiplying sinh(7(fc)) in (4.43) is 
itself a spin operator in the paired subspace. Raising 
(4.43) to the N*"^ power is, then, straightforward: 



(V'(fc))^ 
= cosh(iV7(fc)) 

-sinh(7V7(/c))(cos((5'(fc))T^(fc) - sin((5'(A:))T^(A:)) 

(4.47) 

Since r^(fc)|0) = -|0) and for M even, exp(i7rM) = 1, 
evaluating the matrix element in (4.8) is an easy matter: 



exp[N{K* + K2)] 

X Y[ {cosh(iV7(fc)) + sinh(iV7(fc))cos((5'(fc))} 



0<fc<7r 



(4.48) 



Checking the behavior of 6'{Tr), we see that this equation 
can be written as 



Thus 



Yl {cos{Nj{k)) + smh{Nj{k)) cos((5'(fc))}^/2 



(4.49) 



J ln(cosh(iV7(fc)) + sinh(iV7(fc))cos((5'(fc))) 



fceOjf(+) 



(4.50) 



This is in fact correct whether M is odd or even. The 
result for , the case in which antiperiodicity is enforced 



in each M-spin row, is a simple variation on (4.50). One 



simply replaces the requirement k e f2jv/(+) in the sum 
by fc e flMi-)- 



vanish, which turns out to be the case on detailed cal- 
culation. The summand in (1.45) becomes Nj, and the 
sum may be implemented as a contour integral, following 
the methods of ESI 121: 



lim TV^MnZ 

7V-foo 



M 
An 



1 



Ml 



Moj) (4.51) 



where C is the rectangle in the complex w plane 
with sides which are segments of the lines Re(a;) = 
±7r, Im(w) = ±e. The number e is chosen small enough so 
that the only singularities of the integrand inside the con- 
tour are the zeros of the denominator. Equation (1.46) 
can then be written as: 



lim N-^nZf 



M 



duj 7(w) 



dijj- 



1 + e 



(4.52) 



The first term on the right hand side is the bulk term. 
The Casimir contribution is contained in the second term 
on the right hand side of (4.52). The function 7 has 



branch cuts on the imaginary axis symmetrically posi- 
tioned about the real axis. Using the 27r periodicity of 



the integrand in (4.52), the contour in (4.52) will be de- 



formed into one surrounding the branch cut in the upper 
half plane. It is convenient to make the transformation 
uj — ij{u) and interchange Kl and K2 in Eq. (4.44). 
Expanding the denominator, integrating by parts and re- 
summing gives the Casimir free energy from (4.52) as 



F,, 



knT 



du In 



1 + exp 



(^-Mj(u)j (4.53) 



The resulting Casimir force is denoted fcasicyl-) and is 
just: 



fcas{cyl. 



keT 
Att 



du 7(u) [1 - tanh (M7(u)/2)] 



(4.54) 

Note that this expression is invariant under interchange 
of Kl and K2; that is, there is a dual symmetry; further, 
this force is attractive. The scaling form (with Ki = K2 
) is accessed by noting that ^{u) ~ [7(0)^ -I- u^]^/"^ for 
small u and 7(0). We then introduce scaling variables 
y — Mu/2 and x — M^{0)/2 and take the scaling limit, 
giving: 



B. Periodic and antiperiodic boundary conditions 

The limit as iV — > 00 is easily taken, since M < 00 
implies that the sum is finite. Another way of examin- 
ing this is to go back to (4.12). As — ?> 00, we expect 
only the contribution of the maximum eigenvector to sur- 
vive, this provided its scalar product with |0) does not 



M^f,as{cyl.) 



kuT 



dy{x^+y^f'^ 



1 - tanh ( {x^ + y 



2N1/2 



(4.55) 



This implie s the same scaling function as in (3.15). The 
integral in (4.551 has the value 7r^/12 at a; = 0. 
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The evaluation of is also quite straightforward with 
this procedure; the analogue of (4.52) is: 



lim TV^^Z? 



M 



dio 7(w) 



M 
'2^ 



duj- 



iMt 



(4.56) 



with an additional term sgn(rc — r)7(0), which will not 
appear in the Casimir term, since it is independent of 
M. Note that the only other difference is the sign in the 
denominator of the second integrand; this comes about 
because the sum leading to this term is over periodic wave 
numbers. The Casimir free energy in this case is denoted 
by Kas{<^V^-)'^ it is given by: 



knT 



du In [1 — exp {—M^{u))] 



(4.57) 

Note that this expression is also invariant under inter- 



change of K1 and K2- The analogue of (4.55) for the 
Casimir force with the line of reversed bonds is: 



dy (x' + y') 



2U/2 



coth((x2 +?/2)i/2^ _ 1 (4.58) 



The integral for x = has the value 7r^/6. 



C. Free and fixed boundary conditions 

To calculate the Casimir force for free and fixed bound- 
ary conditions we consider the other limits, i.e. M — ^ oo 
with N finite. This yields 



lim 



1 



J1/-J-00 M 



InZf 



1 

47r 



27T 



In 



cosh(iV7(w)) 



+ sinh(iV7(w)) cos 6'{uj) 



duj 
(4.59) 



Extracting the bulk term we are left with an incremental 
part, denoted fy"". This remainder is given by 



peas ^ keT 
f 



In 



(l + cos(5'(w))) 



+ {1 - cos((5'(a;))) cxp[-2A^7(cj)] 



duj 



(4.60) 



Taking the derivative of the above result with respect to 



rpcas 

- 



7(0;) duj 



1 



l + [(l + cos((5'(a;)))/(l 



cos(,5'(a;)))]exp[2Ar7(w)]' 
(4.61) 



and introducing scaled variables as in Section [III B| with 
M replaced by N , we obtain the following expression for 
the Casimir force 



free 
cas 



/OO 
dn^/f^ + fi2 
-00 



(1 



(1 



(4.62) 



which agrees with the expression (3.24) derived by other 
methods. 

An analogous argument allows us to evaluate the 
Casimir free energy for the strip with plus spins at both 
ends. The right hand side of (4.15) can be evaluated in 



a convenient form for calculations by relating it to the 
limit of V when — ^ 0. Following this with the limit 
Kq — > 00 then yields 



lim exp 



M 

E( 

m—l 



I 



-1 



(4.63) 



The states on the right hand side of (4.72) are defined 
by lim/f._j.o |*&±) = \^±)- Thus they are eigenvalues of 
V2(±) written in a form which is useful for the calcula- 
tion. Then a straightforward calculation gives 

($°|(V(±))^|$0) 

Yl [cosh(7V7(fc)) -|-sinh(7V7(fc))coS(5*(fc)]^/^ 



feeOM(±) 



(4.64) 



where the angle S* is defined by ( 4.45 ) and ( 4.46 ) with Kl 
and K2 interchanged, which is just what we could have 
anticipated from duality, as explained above. Evaluating 
the analogue of (4.53) we find 



+ 

cas 



(1 



(1 



) + (! + 



(4.65) 
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Thus 

/c+atW=/caT(-i) (4-66) 

which is a manifestation of duahty. 

V. CONCLUSION 

In this paper, we exploit the exact solutions for the 
partition function of the Ising model in one and two di- 
mensions to examine the impact of boundary conditions 
and dimensionality on the Casimir forces. Our results 
are effectively summarized in Figs. [T]and[2j As both fig- 
ures unambiguously illustrate, boundary conditions alone 
do not suffice to predict the relative strengths of critical 
Casimir forces. This is perhaps most strikingly evident 
in the case of the forces associated with periodic and 
free boundary conditions. While free boundary condi- 
tions give rise to a vanishing critical Casimir force in one 
dimension, in the two dimensional case the same bound- 
ary conditions result in a significant attractive force. By 
contrast, periodic boundary conditions yield a non-zero 
critical Casimir force in one dimension-that is, a force 
that is infinitely larger than the one associated with free 
boundary conditions. On the other hand, in two dimen- 
sions, the Casimir force generated by the same boundary 
conditions is smaller in amplitude than the Casimir force 
resulting from free boundary conditions. 

In the case of anti-periodic boundary conditions, the 
critical Casimir force is repulsive in both one and two 



dimensions. An additional contrast between periodic and 
antiperiodic boundary conditions is evident in the limit 
T = in one dimension; at that point, the Casimir force 
is zero for the case of periodic boundary but is finite in 
the anti-periodic case. 

The methods presented in this paper can also be 
extended to calculate Casimir forces with mismatched 
boundary conditions. 
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Appendix A: Derivation of exact partition functions 
for the Ising Model with antiperiodic-antiperiodic 
boundary conditions 

Here, we make use of results in Ref. [IH| for the reduced 
partition function of the two dimensional square Ising 
model with antiperiodic boundary conditions in both di- 
rections (see Eq. (79) in Ref. [18] ) 



1 



-ill 



sgn 



with = ksT/J and 



(Al) 



N-l M-1 

n n 

p=0 q=0 . 



Ao - Al cos — \^ - A2 cos — — - - A3 cos ' ^ ' ^ ' 



N 



M 



N 



M 



1/2 



(A2) 



Given and thus, 

t==tanh/3J (A3) 

for a square lattice with uniform interaction strength J, 
we find 

Ao = [l+tY (A4) 

Al = 2t{l-t^) (A5) 

A2 = 2t(l-t2), (A6) 

A3 = (A7) 
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N-l M-l 

n n 

p=0 q=0 



(1 + _ 2t{l - t') cos + - 2t(l - t') cos '"('^ + 



N 



M 



(A8) 



The quantity Qq,o in Eq. ( A2 1 needs special attention 



Consider the p = q = term in the product for flo,0: 

[(l+t2)2_4t(l_i2)]l/2 ^ [(l_2t-t2)2]l/2 

= |l-2t-t2| (A9) 



Making use of ( A3 ) , we have for the last line of ( A9 1 



|l-2tanh/3J-tan2/3J| 

= \ 1 cosh^ /3 J - sinh^ /3 J - 2 cosh /3 J sinh (3J\ 

cosh f3J 



— ^|1- sinh 2/3 J| 
cosh pj 



According to (3.10), the critical temperature is at 



cosh f3J coth f3J — 2, 
which can be rewritten as 

sinh 2/3 J = 1 



(AlO) 

(All) 
(A12) 



Thus, the p = q ~ term in the product for iloo is 
proportional to |6'— • The term going as the sign of that 
expression guarantees that the product goes as 9— 9c- It is 
a straightforward exercise to verify that all other factors 
in the product are free of singularities in the vicinity of 
the critical point. 



1. Sum over the variable q 



Now we focus on the sum over the variable q in Eq. 



(A8). There are two cases of interest: v — and i> = 1/2. 
We will analyze the case i/ = making use of the method 
outlined in Section |III C[ The generic sum of interest is 



M-l 

E 

q=0 



6cos?^ 



(A13) 



which is equivalent to the integral 



1 



Mz 



Af-l 



-dz 



J a-b{z + l/z)/2 ~zM - 1 
1 2 r 1 Mz^ 

2THb J z2 - 2az/6+ 1 z*^ - 

1 2 r 

'2wibJ [z - a/b - ^/{a/bY - l){z - a/b + VW^)^ - 1) 



-dz 



1 



1 



- 1 



-dz 



(AM) 



The integral is around the contours shown in Fig. |4] This 
integral is evaluated by closing it around the two poles 
in the first term in the intregrand on the third line of 
(A14). We find after evaluating residues that the result 
of this integration is 




(A15) 



Setting a/b + \J {ajVf- — 1 = F and using the fact that 

we can write 



alb- ^{a/bf - 1 = {a/b+ ^W^F 
( A15 1 as follows 



1)- 



M F^^ + 1 



M 



~ 1 

pAf/2 _|_ p-M/2 



M 



y2 rA^/2 - F 

M 



- 62 



coth 



-M/2 



InF 



(A16) 



14 



The sum of interest is in fact, 

M-l , 

^ In ia-h 

g=0 ^ 



cos — — 
M ) 



(A17) 



which can be obtained by integrating the sum in (A13) 



with respect to a. Making use of the fact that 
|^ln(a/6+V(«/&)2-l) =l/v^ 



T? — }p- , we find the an- 



tiderivative of the last hne of (A16) with respect to a 

sinh I ylnr j (A18) 



2 In 



Finally, if we define 7^ = In Fj, we have a result that be- 
gins to look like a contribution to the expressions utilized 
by Ferdinand and Fisher in the case of periodic boundary 
conditions. To complete the connection, we note that in 
this case 

a = (l + t2)2_2t(l-i2)^og27r(p + Ai) ^^^^^ 



TV 



b = 2t{l-t^) 



(A20) 



Then, 

r„ 



(i + t 



2\2 



i(l-<2) 



27r(p + ^i) 
N 



cosh 2/3 J coth 2/3 J — cos 



N 
N 



cosh 2/3 J coth 2PJ ~ cos ^^^^ | _ 1 



(A21) 



Eq. (A18) corresponds to the logarithm of the product 
over q of f2^o • Reconstituting the product over p, we have 



fj,0 



n ^^"^^ ( 



p=0 



(A22) 



Let's focus on one term in the product for floo: the p = 
term. It is 



sinh ( Y In ( cosh 2^J coth 2/3 J - 1- 



(cosh2/3Jcoth2^J - 1) - 1 



(A23) 



The function cosh2/:JJcoth2/JJ has a minimum at /3J — 
cosh"^ ( W^i (1 + 72) y We write 



(A24) 



/3J - cosh-i (1 + ^)) 



Then, the argument of the hyperbolic sine in ( A23 ) be- 
comes 



M 



In [1 + 4t -t- Sr^ -|- O (r^)] = 2Mt 



(A25) 



which means that as the temperature passes through the 
critical temperature, the hyperbolic sine function changes 
sign. This is the source of the "sgn" term in (3.16) and 

The calculations utilized to obtain fi^g can be adapted 
to the evaluation of the sum over o in f2,,i. After a set 
of tedious but straightforward calculations, we obtain for 
that result 



M 2 



N-1 

JJ^ cosh 

p=0 



M 



(A26) 



Appendix B: Contributions of the terms in (3.161 to 



the partition function in the case of antiperiodic 
boundary conditions 

We first focus on the sums involving hyperbolic cosines. 
We have 



^ Incosh -Af72r+i = ^ / Incosh ^7(2) 



Nz 



N-l 



-dz 



and 



(Bl) 



N-l ^ 

In cosh - Mj2r 

r=l 



1 

27h 



M , Nz^-\ 

In cosh — 7 2; — dz 

2 z^ -1 

(B2) 



The contours of integration in (B2) shown in Fig. |4j the 
contours of integration in (Bl| are identical except for 



a shift along the unit circle. The contours enclose the 
roots of the function l/{z^ — 1), as depicted on the left 
hand side of that figure. They are then deformed to lie 
just inside and outside the unit circle — see the right hand 
side of Fig. |4] We have already established in Section 
III A that the other singularities in the integrands in (Bl ) 



and (B2) lie along the real axis. They are branch cuts 



which are close to the unit circle in the vicinity of the 
critical temperature. 

We now concentrate on the asymptotic expression in 
the vicinity of the critical point and perform expansions 
like those in Section III A and in Eq. (A25l and find the 
following expression 



M 



In cosh ^ V + 



(B3) 



with (f)p — 2Tr{p + 1/2) /N in case of Eq. (Bl I and with 



Pp = 2Trp/N for Eq. p32| ). 

This holds when the complex variable z is in the im- 
mediate vicinity of z = 1, where the unit circle passes 
through the real axis. In fact, (j> = zLi{z— 1). The branch 
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FIG . 4: Contours for the integral on the right hand side of 
(B2l. 




points on the real axis are associated with the values of 
z for which the argument of the natural logarithm goes 
through zero. This occurs when 



cosh ^ \/ 4r2 



A2 - 







(B4) 



We cause the hyperbolic cosine to go through zero by 
making its argument imaginary so that it becomes the 
trigonometric cosine. To do this, we have — > ±i$. 
When (f> > 4t, 



cosh — \l 4r2 



02 



cos — V'l' 



4t2 



The right hand side of (B5| goes to zero when 
M 



2 



Solving for $, we find 



$ = 



-4r2 



(B5) 



(B6) 



(B7) 



Of course, the cosine goes to zero at an infinite 
set of points on the real axis, corresponding to 
(M/2)V$2 - 4t2 = (2n + l)7r/2. That is to say, there 
is an infinite number of branch points on the real axis. 
However, we are interested in the one closest to the unit 



circle, corresponding to the conditions in (B6) and (B7|. 
The two branch points of interest are just inside and just 
outsize the unit circle. The contribution of the integral 
around the branch cut ending at those branch points will 
be dominated by a multiplicative exponential factor. The 
portion of the distorted contour that lies just inside the 
unit circle can be pushed towards the origin so that it 
wraps around the branch cut. The general strategy for 
deforming the contours is illustrated in Fig. [5j Note that 
the branch cuts are along the real axis. To assess the con- 
tribution of the contour lying just outside the unit circle, 
we rewrite the last portions of the integrands as follows 



(B8) 
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FIG. 5: The three contours over which integrations are per- 
formed. The dashed circle is the unit circle, and the circular 
contour actually lies right on the unit circle. 



Therefore, the right hand side of (Bl) becomes 



/In cosh — 7(z) — dz 
2 z 



(BIO) 
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and that of (B7) becomes 
M 
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1 



1 



dz (B12) 



The first terms in (Bll) and (B12) are integrated around 
the unit circle and become the dominant contribution 
to the sum we are performing. The second terms are 
deformed to wrap around the branch cut that is outside 
the unit circle. 

For the branch cut just outside the unit circle, z = 
1 -I- $, and the factor of interest goes as 



= (1 + $)-^ 



(B13) 



In the above, we have neglected the difference between 
factors going as (1±<I>) and 1. This is because $ is going 
to be smaU (0(1/M)). 

Looking at the contour around the branch cut outside 
the unit circle, we see that the magnitude and the sign 
of the contribution is fundamentally the same for both of 
the In cosh sums. On the other hand, because of the sign 
differences in the second terms on the right hand sides 
of ( B8 1 and ( B9 ) , we see that the two sums will differ by 



exponentially small amounts, the size of the difference 
being dominated by Q-^^/^^^TMYA^ir^ ^ ^.j^g ^^^^ ^-j^g 
branch cut just inside the unit circle, z = 1 — $, and we 
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have to consider the factor 



N 



= (1-$) 
^ gArin(l-*) 



(B14) 



Exponentiating the sums and expanding the exponen- 
tials with respect to the very small (in fact infinitesi- 
mally small in the thermodynamic limit N — ^ oo) terms, 
we find that there is a remainder after the subtraction 
that goes as 



exp 



TV V4r2 + (7r/M)2 



/ (In cosh — y^T^^T^ 



(B15) 



The second term in the exponential comes from the inte- 
gral around the unit circle. There are also multiplicative 
terms, and they may depend non-trivially on key vari- 
ables. However, they have only marginal effect on the 
overall amplitude of the expression in (B15). 



In the case of the hyperbolic sine sums, the analysis 
is basically the same, up to the point at which we look 
for the branch point. Here, we look for the values of 
z at which the hyperbolic sine goes through zero. This 
happens when 



sin ^ _ 4^2 ^ 



(B16) 



The value of <& for which this occurs that is closest to 
the intersection between the unit circle and the real axis 
{z = 1) is 

$ = ±2r (B17) 

The same set of steps taken above leads us to a remainder 
if the exponentiated hyperbolic sine sums are subtracted 
from each other that goes as 



exp 



-N\A7 



N 

271 



M I 

In sinh — \ 
2 ^ 



The term in (B18) dominates the term in (B17). To 



d<j) 

(B18) 



see this, let's focus on the partition function in the region 
of interest, writing 



t 

M 
M 



(B19) 
(B20) 



Then, if we take the difference between the exponents in 
(IB 171) and (|B18|), we have 



iV 
M 



- / Incoth y^t^+n^ dn 



V^2 + 4f2 + 2|i| 



(B21) 



The plot of this difference is shown in Fig. [6] Note that 




FIG. 6: The difference between the exponents in (B17l and 
( |B18| |, plotted versus the scaled reduced temperature, t. 



it is always negative. 



Appendix C: Duality and boundary conditions 

Establishment of the duality between the two dimen- 
sional Ising lattice with free boundary conditions and the 
lattice with fixed boundary conditions can be achieved in 
two ways: by comparison of high and low temperature 
expansions for the partition function or with the use of 
automorphisms in the Kauffman operators. We will re- 
view both approaches in this Appendix 



1. Duality via high and low temperature 
expansions 

The first demonstration of Kramers- Wannier duality 
relation between the Casimir force for fixed spin and free 
boundary conditions [351 US] utilizes a standard argu- 
ment based on high and low temperature expansions [2 7) . 
Figure [T] shows the Ising lattce subject to fixed bound- 
ary conditions along with its dual. The fixed-boundary- 
condition lattice has bonds represented as solid black 
lines, and the sites at which the spins sit are shown as 
black circles. The expansion that we consider for the par- 
tition function of this lattice is the low temperature one, 
in which we assume a "ground state" in which all spins 
are up, and we expand with respect to down spins. This 
generates a graphical expansion for the partition function 
containing terms such as the one displayed in Fig. |8] 

Figure [9] shows one of the sections of reversed spin in [8] 
The red dashed lines denote the boundary of that region. 
As indicated in the figure, the strength of the horizontal 
bonds is Ki and the strength of the vertical bonds is K2 ■ 
The contribution to the partition function of a portion 
of the boundary intersecting a horizontal bond is e~^^^ 
corresponding to the alteration of the Boltzmann factor 
associated with the mismatch of a pair of spins; when 
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FIG. 7: The lattice with fixed spin boundary conditions 
(black) and the dual lattice with free boundary conditions 
(red). 



that all contributions to the low temperature expansion 
consist of complete closed boundaries 

The next step in the development is to generate a high 
temperature expansion for the partition function of the 
spin system on the dual lattice. The spins on this lattice 
interact via nearest neighbor bonds that are both hori- 
zontal and vertical as in the original lattice. The strength 
of the horizontal bonds is and of the vertical bonds 
is The high temperature expansion is based on the 
decomposition 



exp[Kaa'] = cosh K + aa' sinh K 

= coshKil + aa'tanhK) (CI) 

The expansion is in the second term in parentheses on 
the last line of (CI ). Because there is precisely one such 



term for each pair of spins and because each spin must 
appear an even number of times in the expansion, the 
high temperature expansion is represented graphically by 
figures like the border shown in Fig. [9j this time with the 
portions of the figure standing for bonds between spins 
on the dual lattice, as shown in Fig. [Toj The portion is 
shown in Fig. flO) Each segment of the portion yields a 



FIG. 8: Graphical representation of a term in the low tem- 
perature expansion of the partition function. The spins on 
sites inside the regions enclosed by dashed red lines point in 
a direction opposite to that in which all other spins point. 
The spins on the boundary, which are fixed, are not allowed 
to reverse. 



K- 



Ko . 




FIG. 9: A region of reversed spins 

the portion of the boundary intersects a vertical bond 
the contribution is e~^^^ . 

An important feature of this expansion is that fixed 
boundary conditions, in which the spins at the periph- 
ery of the system are constrained to point up, guarantee 



FIG. 10: The portion of the high temperature expansion cor- 
responding to the detail in the low temperature expansion 
shown in Fig. [9] 



factor tanhKl if it is horizontal and tanhXj. Perfect 
correspondence is achieved if 

UnhKl = exp[-2ii'2] (C2) 
tanhXa* = cxp[-2Ki] (C3) 

This is the standard set of duality relations. 

Here, in order that the high temperature expansion 
also consist of closed contours, we must allow all spins 
on the dual lattice to vary without restrictions. That is, 
the boundary conditions on the dual lattice must be free. 
This establishes the complete duality between the two 
versions of the two dimensional Ising model. 
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